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Studies of low-frequency resistance noise show that the glassy freezing of the two-dimensional (2D) 
electron system in the vicinity of the metal-insulator transition occurs in all Si inversion layers. The 
size of the metallic glass phase, which separates the 2D metal and the (glassy) insulator, depends 
strongly on disorder, becoming extremely small in high-mobility samples. The behavior of the second 
spectrum, an important fourth-order noise statistic, indicates the presence of long-range correlations 
between fluctuators in the glassy phase, consistent with the hierarchical picture of glassy dynamics. 
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Despite many theoretical and experimental efforts, the 
metal-insulator transition JlJ (MIT) in two-dimensional 
(2D) systems remains controversial. Since the apparent 
MIT occurs in the regime where both electron-electron 
interactions and disorder are strong, it has been sug- 
gested that the 2D system undergoes glassy ordering in 
the vicinity of the MIT. The proposals include freezing 
into a Coulomb ||, ||, Wigner or spin glass ||. In- 
deed, a recent study of low-frequency resistance noise in 
an extremely low-mobility (high disorder) 2D electron 
system in Si demonstrated |}| glassy freezing, which oc- 
curred in the metallic phase as a precursor to the MIT. 

Here we report a detailed study of low-frequency re- 
sistance noise in a 2D electron system (2DES) in Si in 
the opposite limit of very low disorder, where the metal- 
lic drop of resistivity p with decreasing temperature T 
is most pronounced. Such samples have been studied 
extensively |j} Q in the context of a 2D MIT using mag- 
netotransport measurements. We find that, similar to 
the case of low-mobility samples, the behavior of several 
spectral characteristics of noise in these high-mobility de- 
vices indicates a sudden and dramatic slowing down of 
the electron dynamics at a well-defined electron density 



metallic glass phase (n c < n s < n g ), which separates the 
2D metal and the (glassy) insulator, is consistent with 
recent predictions of the model of interacting electrons 
near a disorder-driven MI T J9| . Furthermore, by analyz- 
ing the second spectrum |l0[ pd|, an important fourth- 
order noise statistic, we have established the presence of 
long-range correlations between fluctuators in the glassy 
phase, which provides an unambiguous evidence for the 
onset of glassy dynamics at n g . The results are consis- 
tent with the picture in which noise, in the glassy phase, 
results from transitions between many metastable states 



with a hierarchical structure |12 with each transition be- 
ing a reconfiguration of a large number of electrons. 

Measurements were carried out on n-channel Si MOS- 
FETs with the peak mobility p w 2.5 m 2 /Vs at 4.2 K, 
fabricated in a Hall bar geometry with Al gates, and oxide 
thickness d ox = 147 nm [Q. The resistance R was mea- 
sured down to T = 0.24 K using a standard four-probe 
ac technique (typically 2.7 Hz) in the Ohmic regime. A 
precision DC voltage standard (EDC MV116J) was used 
to apply the gate voltage, which controls n s . Contact re- 
sistances and their influence on noise measurements were 
minimized by using a split-gate geometry, which allows 



corresponding to the transition to a glassy one to maintain high n s (w 10 cm ) in the contact 



phase. Since the two sets of devices, Si metal-oxide- 
semiconductor field-effect transistors (MOSFETs), differ 
considerably by their peak mobility, which is a rough 
measure of the disorder, and span essentially the entire 
range of Si technology, we conclude that the observed 
glass transition is a universal phenomenon in Si inver- 
sion layers. The experiments, however, have also revealed 
an important difference between low- and high-mobility 
samples. In low-mobility devices, n g ~ 1.5 n c ||, where 
n c is the critical density for the MIT determined from 
the vanishing of activation energy and the tempera- 
ture coefficient of p changes sign at n* > n g > n c . In 
high-mobility structures, on the other hand, the onset 
of glassy dynamics seems almost to coincide with the 
MIT, i. e. n g n c w n*. It is interesting that the ob- 
served strong dependence on disorder of the size of the 



region while allowing an independent control of n s of the 
2D system under investigation in the central part of the 
sample (120 x 50 pm 2 ) (Fig. [I] inset). Nevertheless, care 
was taken to ensure that the observed noise did not come 
from cither the current contacts or the regions of gaps in 
the gate. For example, since the noise measured across a 
resistor connected in series with the sample and having a 
similar resistance was at least three times lower than the 
noise from the central part of the sample, the effect of the 
contact noise on the excitation current I exc could be eas- 
ily ruled out. Similarly, the resistance and the noise mea- 
sured between the voltage contact in the region of high 
n s (e. g. #5 in Fig. [l] inset) and the one in the central 
part (#6) were much smaller than those measured be- 
tween contacts in the central part (e. g. #6 and #7). In 
fact, they were in agreement with what is expected based 
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FIG. 1: (p) vs. T for n a (10 10 cm" 2 ) shown on the plot. 
Insets: sample schematic (dashed lines represent gaps in the 
gate), and activation energies us. n s ; n c ~ n*. 



on the geometry of the sample, which proves that the 
submicron gap regions did not contribute to either the 
measured resistance or noise. In order to minimize the 
influence of fluctuations of both I exc and T, some of the 
noise measurements were carried out with a bridge con- 
figuration [[[3) . The difference voltage was detected using 
two PAR124A lock-in amplifiers, and a cross-spectrum 
measurement was performed with an HP35665A spec- 
trum analyzer in order to reduce the background noise 
even further . The output filters of the lock- in ampli- 
fiers and/or spectrum analyzer served as an antialiasing 
device. Most of the noise spectra were obtained in the 
/ = (10 -4 -10 -1 ) Hz bandwidth, where the upper bound 
was set by the low frequency of I exc , limited by the low 
cut-off frequency of RC filters used to reduce external 
electromagnetic noise as well as by high R of the sample. 

Fig. [I] shows the time-averaged resistivity (p) as a 
function of T for different n s ; d(p)/dT = at n* w 
9.7 x 10 10 cm~ 2 . For the lowest n s and T, the data 
are described by the simply activated form (p) oc 
exp(EA/ksT). The vanishing of Ea is often used as a 
criterion to determine n c @. This method (Fig. | inset) 
yields n c ss n* , in agreement with other studies £ ] . 

Figures ||(a) and (b) show the time series of the relative 
changes in resistance AR(t)/(R), and the corresponding 
power spectra Sn(f) oc 1// Q , respectively. In order to 
compare the noise magnitudes under different conditions, 
the spectra were averaged over two octaves [(0.5 — 2) x 
1CP 3 Hz] around / = 10~ 3 Hz. The resulting fraction of 
power Sfi(f = 1 mHz) is taken as the measure of noise, 
and its dependence on n s is shown in Fig. |](a) for T = 
0.24 K. The exponent a is displayed in Fig. ||(b), while 
Fig. ||(c) shows the dependence of Sr(/ = 1 mHz) on 
T for several n s . Below T w 3 K, the noise increases 
with decreasing T, but at high n s where d(p)/dT > 0, 
Sr depends rather weakly on both T and n s . In the 
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FIG. 2: (a) AR/(R), and (b) the corresponding power spectra 
Sn(f), for n s (10 10 cm~ 2 ) shown on the plots; T = 0.24 K. In 
(a) traces are shifted for clarity. In (b) Sn(f) are averaged 
over octaves and multiplied by /, so that 1/f spectrum is 
horizontal on this scale. Solid lines are linear least-squares 
fits. 
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FIG. 3: S R (f) vs. n s (10 10 cm- 2 ) (a), and vs. T (c). (b) a 
vs. n s . Sn(f) has been corrected for the white background 
noise. The critical density n c is shown in (a) and (b). 



vicinity of n c , however, a dramatic change in the behavior 
of Sr is observed. The noise amplitude starts to increase 
strongly with decreasing n s , and a rises rapidly from « 1 
to w 1.8 [^5). This shift of the spectral weight towards 
lower / indicates a sudden and dramatic slowing down of 
the electron dynamics at n g w 10 x 10 10 cm -2 , which is 
attributed to the freezing of the electron glass. The same 
qualitative behavior, together with other manifestations 
of glassiness (slow relaxations, history dependence), was 
observed || in Si MOSFETs with a much higher amount 
of disorder: p was a factor of 40 lower than in the samples 
studied here, and not surprisingly, n g was almost an order 
of magnitude higher. We note that the two sets of devices 
also differ substantially by their geometry, size, and many 
fabrication details (see Refs. pi |?]), all of which leads us to 
conclude that the observed glass transition is a universal 
phenomenon in Si inversion layers, at least in those with 
conventional (d{p)/dT > 0) metallic behavior |Tg| . 

In addition to affecting the values of n c , n g , and n*, 
the disorder clearly plays another, nontrivial role. In par- 
ticular, n c and n g were found to differ from each other 
considerably in low- mobility devices (n c , n g , and n* were 
5.0, 7.5, and 12.9, respectively, in units of 10 n cm~ 2 ) ||, 
whereas in high-mobility devices n g is at most a few per- 
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cent higher than n c [see Figs. ||(a), (b)]. Therefore, the 
emergence of glassy dynamics here seems almost to coin- 
cide with the MIT. Obviously, the size of the intermediate 
(n c < n s < n g ) glass phase depends strongly on disorder, 
in agreement with theoretical predictions Q . 

Earlier studies of noise in R in Si MOSFETs were per- 
formed mostly at T > 4.2 K and n s > 10 12 c m - 2 0, fjjj. 
The observed random telegraph and 1/f noise were at- 
tributed to charging and discharging of traps in the ox- 
ide close to the Si/Si0 2 interface, leading to dS^/dT > 
j fL9| . On the other hand, studies carried out at lower 
n s and T demonstrated |l| clearly that the observed 
1/f noise was an intrinsic property of the conduction in 
a 2D channel and not due to charge trapping. Moreover, 
dSnJdT < was found || for T = L5 > 4 - 2 K - 0ur mea- 
surements at much lower T reveal a dramatic increase 
of Sr with decreasing T . This rules out models of ther- 
mally activated charge trapping jll, 19, 2^| , noise gener- 
ated by fluctuations of T 123] , and a model of noise near 
the Anderson transition |24], as possible explanations. 
Likewise, the models of noise in the Mott and Efros- 
Shklovskii variable-range hopping regimes |25| do not de- 
scribe the data because they predict either dSn/dT > 
or a saturation of Sr below 10-100 Hz, both in clear 
disagreement with the experiment. Therefore, the ob- 
served noise cannot be a result of single electron hops 
even when Coulomb interactions are included through 
the Coulomb gap. We note that no such low-frequency 
saturation of Sr was found in computer simulations of a 
Coulomb glass, where 1/f noise was a result of transitions 
between many metastable states, with each transition be- 
ing a reconfiguration of a large number of electrons [ p6[ . 

We have established that the exponent a ~ 1 in the 2D 
metallic phase (above n g ) in both low- and high-mobility 
samples. On the other hand, a w 1.8 in the glassy phase, 
similar to a found in some spin glasses |27|, p8fl , and sub- 
micron wires in the quantum Hall regime ||29|^ In general, 
such noise with spectra closer to 1/f 2 than to 1// is typi- 
cal of a system far from equilibrium, in which a step does 
not lead to a probable return step. Such high values of 
a may be also obtained if noise results from a superposi- 
tion of a small number of independent two-state systems 
(TSS) 0, |§. However, even though some distinguish- 
able discrete events can be seen at low n s [Fig. ||(a)] , they 
do not show the characteristic repetitive form of stable 
TSS. On the contrary, both the shape and the magni- 
tude of noise exhibit random, nonmonotonic (which ex- 
clude aging) changes with time. A quantitative measure 
of such spectral wandering is the so-called second spec- 
trum 52(/2, /), which is the power spectrum of the fluc- 
tuations of SR(f) with time JToj ] , i. e. the Fourier trans- 
form of the autocorrelation function of the time series of 
SR(f). If the fluctuators (e. g. TSS) are not correlated, 
•52 (/2,/) is white (independent of / 2 ) |h], 11 and equal 
to the square of the first spectrum. Such noise is called 
Gaussian. On the other hand, S2 has a nonwhite char- 
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FIG. 4: (a) Second spectral density S2(/2) fi for 

n s (10 10 cm -2 ) shown on the plot; ft = 1 mHz. In- 
set: log[S2(/2)] for spin glasses Cuo.91Mno.09 [0 and 
Cdo.93Mno.07Te P7|. Solid lines are fits. Exponent 1 — /3 vs. 
n s (10 10 cm -2 ) for (b) high-mobility and (c) low-mobility sam- 
ples. The low-mobility device is the same as the one studied 
in Ref. ^| 



acter, S% oc l/f%~^, for interacting fluctuators |l0|, [TT| . 
Therefore, the deviations from Gaussianity provide a di- 
rect probe of correlations between fluctuators. 

We investigate S2 using digital filtering Q in a given 
frequency range / = (/z,,/h) (usually f H = 2/ L ). 
The normalized second spectra, with the Gaussian back- 
ground subtracted, are shown in Fig. ||(a) for two n s , 
just above and just below n g . It is clear that there is 
a striking difference in the character of the two spec- 
tra. Similar differences are observed between various spin 
glasses (Fig. [I](a) inset), where 62 is white [^7| in the ab- 
sence of long range interactions, and nonwhite [ [To] when 
long range RKKY interaction leads to hierarchical glassy 
dynamics |t^. A detailed dependence of the exponent 
(1 — j3) on n s has been determined for both high- and 
low-mobility samples (Figs. ||(b) and (c), respectively). 
In both cases, S2 is white for n s > n g , indicating that 
the observed 1/f noise results from uncorrelated fluctua- 
tors. It is quite remarkable that S2 changes its character 
in a dramatic way exactly at n g in both types of samples. 
For n s < n g , S2 is strongly non-Gaussian, which demon- 
strates that the fluctuators are strongly correlated. This, 
of course, rules out independent TSS (such as traps) as 
possible sources of noise when n s < n g . In fact, a sud- 
den change in the nature of the fluctuators (i. e. cor- 
related vs. uncorrelated) as a function of n s rules out 
any traps, defects, or a highly unlikely scenario that the 
observed glassiness may be due to some other time depen- 
dent changes of the disorder potential itself. Instead, it 
provides an unambiguous evidence for the onset of glassy 
dynamics in a 2D electron system at n g . 

In the studies of spin glasses, the scaling of £2 with re- 
spect to / and /2 has been used |l(| to unravel the glassy 
dynamics and, in particular, to distinguish generalized 
models of interacting droplets or clusters (i. e. TSS) 
from hierarchical pictures. In the former case, the low-/ 
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FIG. 5: Scaling of & measured at frequency octaves / = 
(/i,2/_L) for (a) high-mobility and (b) low-mobility samples. 



noise comes from a smaller number of large elements be- 
cause they are slower, while the higher-/ noise comes 
from a larger number of smaller elements, which are 
faster. In this picture, which assumes compact droplets 
and short-range interactions between them, big elements 
are more likely to interact than small ones and, hence, 
non-Gaussian effects and w iU be stronger for lower /. 
^(/b, /), however, need to be compared for fixed fz/f, 
i. e. on time scales determined by the time scales of the 
fluctuations being measured, since spectra taken over a 
fixed time interval average the high-frequency data more 
than the low-frequency data. Therefore, in the inter- 
acting "droplet" model, S^fo, /) should be a decreas- 
ing function of / at constant fil f . In the hierarchical 
picture, on the other hand, 5*2 (/2, /) should be scale in- 
variant: it should depend only on f2/f, not on the scale 
/ [0 . Fig. H shows that no systematic dependence of $2 
on / is seen in our samples, which signals that the sys- 
tem wanders collectively between many metastable states 
related by a kinetic hierarchy. Metastable states corre- 
spond to the local minima or "valleys" in the free energy 
landscape, separated by barriers with a wide, hierarchical 
distribution of heights and, thus, relaxation times. Inter- 
vallcy transitions, which are reconfigurations of a large 
number of electrons, thus lead to the observed strong, 
correlated, l//-type noise, remarkably similar to what 
was observed in spin glasses with a long-range correlation 
of spin configuration We note that, unlike droplet 
models |3l| , hierarchical pictures of glassy dynamics j3S| 
do allow for the existence of a finite T (or finite Fermi en- 
ergy) glass transition in presence of a symmetry-breaking 
field, such as the random potential in an electron glass. 

In summary, by studying the statistics of low-/ resis- 
tance noise, we have established that the glassy ordering 
of a 2DES near the MIT occurs in all Si inversion lay- 
ers. The size of the metallic glass phase, which separates 
the 2D metal and the glassy insulator, depends strongly 
on disorder, becoming extremely small in high-mobility 
samples. The properties of the entire glass phase are con- 
sistent with the hierarchical picture of glassy dynamics, 



similar to spin glasses with long-range correlations. 
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